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We investigate the thermodynamic properties of Heisenberg ferrimagnetic mixed-spin chains 
both numerically and analytically with particular emphasis on the combination of ferromagnetic 
and antiferromagnetic features. Employing a new density-matrix renormalization-group technique 
as well as a quantum Monte Carlo method, we reveal the overall thermal behavior: At very low 
temperatures, the specific heat and the magnetic susceptibility times temperature behave like T^l'^ 
and r~^, respectively, whereas at intermediate temperatures, they exhibit a Schottky-like peak and 
a minimum, respectively. Developing the modified spin-wave theory, we complement the numerical 
findings and give a precise estimate of the low-temperature behavior. 



I. INTRODUCTION 

Low-dimensional quantum magnets with two kinds 
of antiferromagnetically exchange-coupled centers have 
been attracting much current interest. Several authors 
constructed integrable Hamiltonians and extracted 
suggestive critical phenomena from them. Alternating- 
spin antiferromagnets with singlet ground states stimu- 
late us to study again the nontrivial gap problem Q . Em- 
ploying the nonlinear-cr-model technique, various mixed- 
spin systems such as linear chains |^,^ and ladders |^] 
have systematically been studied with particular empha- 
sis on the competition between massive and massless 
phases. 

In this context a remarkable attention has been di- 
rected to ferrimagnetic mixed-spin chains ||9f-[l^, which 
are the subject in the present article. Since we may ex- 
pect gapless excitations from magnetic ground states, it 
is less interesting there whether the system is massless or 
massive. Using a field-theoretical argument, Alcaraz and 
Malvezzi |^ predicted that mixed-spin isotropic Heisen- 
berg ferrimagnets should exhibit quadratic dispersion re- 
lations. Their prediction was numerically verified and 
the quadratic dispersion was explicitly visualized p4[ . 
Thus quantum ferrimagnets are expected to behave like 
ferromagnets at low temperatures. On the other hand, 
conventional spin- wave calculations ||l^,|l^ and a pertur- 
bation approach p^ ] from the decoupled-dimer limit sug- 
gest that quantum ferrimagnets should exhibit nontrivial 
gapped excitations as well, which have the effect of en- 
hancing the ground-state magnetization and are there- 



fore of antiferromagnetic nature. A quantum Monte 
Carlo (QMC) technique and an exact-diagonalization 
method actually observed the two distinct low-lying 
excitations and showed that the mixed nature remains 
unchanged as long as the model is isotropic. 

Motivated by the revealed low-energy structure, two 
of the present authors have investigated the ther- 
modynamic properties of Heisenberg ferrimagnetic spin 
chains focusing on the idea of coexisting ferromagnetic 
and antiferromagnetic aspects. On the one hand, a QMC 
method was employed to find the overall thermal behav- 
ior. On the other hand, a modified spin-wave (MSW) 
theory was introduced to clarify how the two distinct ex- 
citations contribute to the thermodynamics. The present 
article is a developed richer version of the preceding work. 
In order to inquire further into the low-temperature be- 
havior, we here employ additional numerical tools. A 
quantum transfer-matrix (QTM) method allows us to ob- 
serve how the ferromagnetic character grows with the in- 
crease of the system size. A brand-new density-matrix 
renormalization-group (DMRG) technique helps us to go 
down into the low-temperature region that was never at- 
tained before. Furthermore, taking account of certain 
interactions between spin waves, we refine the MSW the- 
ory so as to bring a more accurate description of the 
thermal quantities at low temperatures. Although it is 
rather hard, even with interacting spin waves, to obtain a 
quantitative description of the overall thermal behavior, 
a grand-canonical approach to ferrimagnets in terms of 
spin waves is suggestive enough and interesting in itself. 
We fully argue how we should modify the conventional 
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spin-wave theory in an attempt to construct the thermo- 
dynamics of quantum ferrimagnets. 

We consider two kinds of spins S and s alternating on 
a ring with antiferromagnetic exchange couphng between 
nearest neighbors, as described by the Hamiltonian 



N 



iV 



4). 



(1) 

where N denotes the number of unit cells, S represents 
a bond alternation, /is is the Bohr magneton, and we 
have set the g factors of both spins S and s equal to 
g. We assume that S > s, which remains general 
enough to describe alternating-spin ferrimagnets. The 
Lieb-Mattis theorem shows that the Hamiltonian (|^), 
unless a field applied, has [2(5 — s)N + l]-fold degen- 
erate ground states. The ferromagnetic and the anti- 
ferromagnetic excitations, which lie in the subspaces of 
M < (S - s)N and M > {S - s)N, respectively, in- 
deed show a quadratic dispersion and a gapped spec- 
trum iQ, where M — X)jLi('S'| + s|) is the total mag- 
netization. The antiferromagnetic gap, namely, the gap 
between the ground state and the lowest excitation to 
M = {S - s)N + 1, was estimated to be 1.75914(1) J in 
the thermodynamic limit. The correlation length of the 
system is so small that it is barely of the length of the 
unit cell [foljl^ even at the Heisenberg point 6=1. 



II. NUMERICAL STUDY 



where n is a Trotter number, /3 = (fceT) ^ with the 
Boltzmann constant fee, and 
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In order to accelerate the thermodynamic calculation and 
refine its accuracy, we make use of the improved algo- 
rithm in Ref. |22|. Calculations are carried out at several 
values of n and N and they are extrapolated to the limit 
of n ^ oo and N oo. Although we have treated chains 
of iV = 24, iV = 32, and N = 48, the size dependence of 
the thermal quantities is so weak that we find no differ- 
ence beyond the numerical uncertainty even between the 
calculations at TV = 24 and = 32 except for very low 
temperatures. Both the specific heat C and the magnetic 
susceptibility x s-re directly evaluated through formulae 

M 



C = kBPmQ2) + (Ql) ~ (Qif) , 
X = <?Vb/3 ((M2) - (Ml)') , 
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with 
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In this section, using several numerical tools, we cal- 
culate the specific heat and the magnetic susceptibility 
at zero field. The spin- wave ||l^,|l^ and the perturba- 
tion ||l^ calculations suggest that the low-temperature 
properties of the model are qualitatively the same regard- 
less of the values of S and s as longas they differ from 
each other. Alcaraz and Malvezzi performed finite- 
size calculations combined with a scaling analysis in the 
cases of (5,s) = (1,1/2) and (5,s) = (3/2,1/2) and in- 
deed concluded that the appearance of quadratic disper- 
sion relations should be expected for arbitrary isotropic 
mixed-spin chains showing ferrimagnetism instead of an- 
tiferromagnetism. Thus we restrict our numerical inves- 
tigations to the case of {S, s) = (1, 1/2). 
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where (A) denotes the thermal average of A at a given 
temperature a. 
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^' and crp'^ are Ising spins on the 



site i at the Zth Trotter layer, and the local Boltzmann 
factors P^2j-P and p'^p are defined in terms of the eigen- 
states of 5? and sf as 



(8) 



A. Procedure 



We employ the QMC method based on the Suzuki- 
Trotter decomposition ||2^ of the checkerboard type . 
The partition function Z = Tr[e~''^] is approximately 
decomposed as 
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L i=l,3, 



(2) 



As far as we replace the trace by an importance sam- 
pling, it is hardly feasible to take grandcanonical aver- 
ages at very low temperatures. In an attempt to avoid 
the difficulty, we may integrate out ^4 25| the (1 -I- 1)- 
dimensional Ising system but then have to give up either 
low temperatures or long chains. We can in principle 
reach an arbitrary low temperature at the expense of the 
system size and such an attempt is to a certain extent 
fruitful in the present system whose correlation length 
is very small. Actually, constructing transfer matrices 
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in the chain direction, we will later show QTM calcu- 
lations of the specific heat. Although the QTM calcu- 
lations themselves are so accurate as to be regarded as 
exact, numerical differentiations of them may introduce 
errors. Here we directly calculate the internal energy, 
rather than the partition function, and numerically differ- 
entiate it once. We keep the final results highly accurate 
by taking raw data at regular intervals of k-oT/ J = 0.01. 

On the other hand, constructing transfer matrices in 
the Trotter direction, we can partly reveal the thermal 
behavior of an infinite chain p6[ |. However, the expo- 
nential growth of the matrix size with the increase of n 
makes our access to low temperatures unfeasible. In or- 
der to overcome the difficulty, we introduce the DMRG 
technique into our investigations. Based on White's orig- 
inal idea developed at T = 0, we apply it to the 
renormalization of transfer matrices ||28| , P9|] instead of 
Hamiltonians. This extension of the DMRG method has 
recently been applied to the thermodynamics of several 
low-dimensional magnets |30-3^] and indeed brought us 
plenty of results at very low temperatures. 

The core idea of the so-called transfer-matrix DMRG 
method can be summarized as reaching large n by re- 
ducing the size of the transfer matrix with the use of the 
DMRG algorithm, which is highly successful in obtaining 
effective low-energy Hamiltonians. We introduce a step- 
width /3o and go successively down to lower temperatures 
— (n/3o)^^ increasing n linearly. The transfer-matrix 
DMRG calculation is technically different from the orig- 
inal procedure in that we have to handle nonsymmetric 
density matrices there. In order to avoid the exponen- 
tial growth of the matrix size, we keep the number of 
states in the density matrix constant at a predetermined 
number m throughout the calculation. At each iteration 
n — > n -|- 1, we first produce the transfer matrix at the 
step n -f 1 by making a matrix product of the transfer 
matrix at the step n and the local Boltzmann factor (H) 
with /3 = n/3o- Next we unite the two transfer matrices 
at the step n + 1 into an augmented matrix of double 
width in the Trotter direction. We then calculate the 
left and right eigenvectors of the enlarged matrix with 
the maximum eigenvalue and construct a density matrix 
from them. It is the eigenstates corresponding to the m 
largest eigenvalues of the density matrix that should be 
retained to give an approximate description of the trans- 
fer matrix at the step n + 1. The thus-obtained transfer 
matrix allows us to calculate the free energy at each tem- 
perature f3~^ = {nf3o)~^. We here carry out the direct 
estimation of the internal energy and the magnetiza- 
tion, which strongly reduces numerical errors in the final 
results, the specific heat and the magnetic susceptibility. 

Obviously the two controllable parameters /?o and m 
decide the precision of the transfer-matrix DMRG calcu- 
lation. The Trotter decomposition is refined as /3o 0, 
whereas the loss of information is reduced as m — > oo. 
At high temperatures, we are compelled to work with 
small Trotter numbers, while few states are discarded. 
At low temperatures. Trotter numbers are large enough, 



whereas many states are lost due to the large number of 
iterations. Thus, the convergence of the calculation pre- 
dominantly depends on f3Q at high temperatures, while 
on m at low temperatures, as was actually observed 

As we have conclusive QMC results at high and in- 
termediate temperatures, we are particularly interested 
in low-temperature findings from the transfer-matrix 
DMRG calculation. We therefore invest computational 
resources mainly in augmenting m. Setting PqJ and m 
to 0.2 and 128, respectively, we have found that at in- 
termediate temperatures the specific heat is somewhat 
overestimated, whereas at low temperatures the preci- 
sion almost reaches three digits. Calculations at m — 80, 
rn = 96 and m — 128 fully converge at T ^ 0.05 and 
almost converge at 0.04 ^ T ^ 0.05. We were not able 
to go further down below successfully. This is mainly 
due to the macroscopically degenerate ground states of 
ferrimagnets, which imply that a huge number of states 
must be retained in our access T ^ 0. 



B. Results 

We show in Fig. ^ the temperature dependence of the 
specific heat of the (5, s) = (1, 1/2) Heisenberg ferrimag- 
netic spin chain. We find in Fig. y(a) a good agreement 
between the QMC and DMRG calculations. At interme- 
diate temperatures, the specific heat exhibits a Schottky- 
like anomaly typical of antiferromagnets. Let us consider 
the model at 5 = in an attempt to clarify the origin of 
this characteristic peak. Now the model is decoupled 
into dimers and the specific heat per unit cell is simply 
obtained from an isolated dimer, which is a pure two- 
level system with doubly degenerate ground states and 
fourfold degenerate excited states separated above by an 
energy gap A = 3J/2. The specific heat is given by 
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r{f3A) 



where r is the ratio of the degeneracy of the excited 
states to that of the ground states. Moving away from 
the decoupled-dimer point, the low- lying states come to 
exhibit dispersion and macroscopic degeneracy. At the 
Heisenberg point S = 1, for N elementary cells, the 
ground state is (A^-|- l)-fold degenerate, whereas the low- 
est antiferromagnetic excited state is {N + 3)-fold de- 
generate We note that the isolated dimer may be 
regarded as the Heisenberg chain of = 1. Now we 
make an attempt to fit the intermediate-temperature be- 
havior to the curve (^) allowing only a single additional 
adjustable parameter A: 
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(10) 



Here we replace r and A by their values at the Heisen- 
berg point, the degeneracy ratio of the lowest antiferro- 
magnetic excited states to the ground states and the gap 
between them, respectively: 
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= 1, /1= 1.759J. 



(11) 



We obtain a fine fit with A = 1.7, which is also shown 
in Fig. 0(a). We stress that the Schottky-type spe- 
cific heat (lOh not only fits the peak but also well de- 
scribes the high-temperature decay. Equation com- 
bined with the condition (11) gives an asymptotic high- 
temperature behavior 1.3(/3J)^, which is well consistent 
with the high-temperature series-expansion result (PJ)'^ 
i pTf . Thus we find that the temperature dependence of 
the specific heat straightforwardly reflects the antiferro- 
magnetic gap and therefore the model can be treated well 
as a two-level system unless it is at low enough temper- 
atures kBT < Z\ ~ 1.759J. 

On the other hand, the DMRG results fully allow us 
to guess the T^/^ asymptotic behavior characteristic of 
ferromagnets. This ferromagnetic feature grows as we 
move toward the Heisenberg point from the decoupled- 
dimcr limit or alternatively as the system size increases. 
From this point of view, we observe the chain-length de- 
pendence of the specific heat in Fig. |l|(b). The naive 
QTM method is most useful in short chains at low tem- 
peratures. We find that the antiferromagnetic feature is 
rather settled even in short chains, whereas the growth 
of the ferromagnetic feature is relatively slow. The ele- 
mentary excitations of ferromagnetic nature are predom- 
inantly related to spin I's, while those of antiferromag- 
netic nature originate from interactions between the two 
kinds of spins . These excitations may be identified 
with the local ones within a unit cell in the vicinity of 
the decoupled-dimer limit, as shown in Fig. |^. Figure 
|l|(b) suggests that the delocalization effect is more es- 
sential to the appearance of the ferromagnetic feature. 
Actually the ferromagnetic excitations constitute a wider 
band than the antiferromagnetic ones |l^ . 

We show in Fig. ^ the temperature dependence of the 
magnetic susceptibility times temperature of the {S, s) — 
(1,1/2) Heisenberg ferrimagnetic spin chain. We find 
again that the QMC and DMRG calculations are in excel- 
lent agreement. The product xT diverges as at low 
temperatures, while it approaches [S'(S'-|-l)-|-s(s-|-l)]/3 — 
11/12 at high temperatures. The low-temperature di- 
vergence is reminiscent of the ferromagnetic susceptibil- 
ity ]35|] . We further note that xT shows a minimum 
while shifting from the quantum ferromagnetic behav- 
ior to the classical paramagnetic behavior. Generally 
xT is a monotonically decreasing function in ferromag- 
nets, while a monotonically increasing function in anti- 
ferromagnets [|3|. The susceptibilities of gapped anti- 
ferromagnets such as Haldane systems vanish exponen- 
tially [^,^ 34 1 . Hence the temperature dependence we 
are observing may be regarded as a ferromagnetic-to- 
antiferromagnetic crossover. Actually the minimum of 
xT appears at a temperature around which the specific 
heat shows the Schottky-like peak. 

In sum, a combination of ferromagnetic and antifer- 
romagnetic features in ferrimagnets has been demon- 



strated. The DMRG technique has allowed us to have 
access to very low temperatures. In an attempt at going 
further down to even lower temperatures, we develop the 
MSW theory for ferrimagnets jl^ in the next section. 
We not only aim at obtaining the precise description of 
the low-temperature behavior but are also interested in 
describing the overall thermal behavior based on the spin- 
wave picture. 



III. MODIFIED SPIN- WAVE APPROACH 



For years the conventional spin-wave theory 
was plagued by the difficulty that the zero-fieU mag- 
netization diverges in low-dimensional magnets. The 
low-tcmperature series expansion within the theory only 
brings us the leading term of the specific heat and noth- 
ing correct for the susceptibility. However, imposing a 
constraint on the magnetization, Takahashi [^|4^ suc- 
ceeded in describing the low-temperature thermodynam- 
ics of low-dimensional ferromagnets. His idea was further 
applied to several quantum antiferromagnets and 
its wide applicability was established. Recently two of 
the present authors |T7[ have introduced this modified 
spin-wave theory into quantum ferrimagnets and demon- 
strated that it is quite useful in understanding their char- 
acteristic features. Here we develop our argument taking 
into account interactions between spin waves. 



A. Dispersion relations 

We have observed that the low-temperature thermody- 
namics well refiects the dispersion relation of the ferro- 
magnetic branch, whereas the intermediate-temperature 
behavior is well attributed to that of the gapped anti- 
ferromagnetic branch. This scenario is indeed supported 
by the MSW theory [|l^. A combination of the ferro- 
magnetic and antiferromagnetic spin waves on which a 
particular constraint is imposed reproduces the T^/^ ini- 
tial behavior, the Schottky-like peak, and the decay 
of the specific heat, and the T"^ divergence and the 
decay of the susceptibility (see Fig. [sjand Fig. ||). How- 
ever, we should be reminded that the conventional spin- 
wave calculation, based on which we have started our 
first attempt JlTf to construct the thermodynamics, con- 
siderably underestimates the antiferromagnetic gap A. 
Consequently, the MSW approach succeeds in reproduc- 
ing the Schottky anomaly itself, to be sure, but fails in 
locating it at a correct temperature. In order to obtain a 
better description, we here refine the original spin-wave 
theory. 

Let us introduce bosonic operators through the 
Holstein-Primakoff transformation fB6l 



S+ - (25-a]a,)i/2aj, 
4 ^b]i2s-blb,y/\ 



(12) 
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where we regard S and s as quantities of the same order, 
that is, 0(5*) = 0(s). The Hamiltonian (§) with 6 = 1 
and if = is expressed in terms of the bosonic operators 
as 



H = Sclass+W0+Wl+O(^"') 



where 



-2sSJN, 



Ho = jYj\'^sa]aj + 2Sb]hj 

j 

+V^^{aj + aj+i)bj + h.c 



(13) 
(14) 

(15) 



+ ^/s/Sa]a'j{bJ + bj-i) + h.c. 



+ (a]aj + a]+iaj+i)6]&„ 



(16) 



The treatment of the quartic interaction ( |16D is not so 
canonical as that of the quadratic Hamiltonian ( p^ ) . A 
variety of approaches are possible |l^. Here, in an at- 
tempt to obtain the dispersion relations beyond the non- 
interacting spin-wave theory, we first diagonalize TLo and 
next extract corrections from TLi. 
The Bogoliubov transformation 



au = coshflfc afe -I- sinh^^ b\ 
Pk = sinh6'fc al + cosh6'fc bk 

combined with 

j 

bi - —\^e-~'^^^+^/^^b- 



(17) 



(18) 



and 



2VSs fk 

cos — 

S + s I 2 



tanh26'fc 
diagonalizes Ho as [ jl0[|12[ 



(19) 



(20) 



where we have taken twice the lattice constant as unity. 
The first term in ( |20| ) is a quantum correction to the 
ground-state energy of order 0(5^), 



k 



(21) 



with 



ujk = \ (S 



4Sssm^{k/2) . 



(22) 



The following terms are the ferromagnetic and antiferro- 
magnetic spin-wave modes of order 0{S^), whose disper- 
sion relations are, respectively, given by 



uJkT{S - s). 



(23) 



The lower-energy mode shows a quadratic dispersion re- 
lation at small fc's, which is consistent with the T^/"^ ini- 
tial behavior of the specific heat. On the other hand, 
the gap between the two branches is exactly J, which 
considerably contradicts the numerical estimate 1.759J. 

Now we pick up relevant contributions to the disper- 
sions, as well as to the ground-state energy, from Tii. 
Employing the Wick theorem, we rewrite Tii as 



■Hi 



El - 



k 



(24) 



where iJinci contains irrelevant terms such as a^Pk and 
Wquart contains residual two-body interactions, both of 
which are neglected in the following. The correction to 
the ground-state energy, Ei , and ones to the dispersions, 
5uj^ , are, respectively, given by 



El = -2JN 



Ft 



Slj^ = 2(5* 4- s) A 



sin^ 



s + 

^2 



'Ss 



ujk±{S - s) 



(25) 
(26) 



with 



^1 = ^E(co^h20''^ - 1) . 



y^cosl — 1 sinh26'i. 

2N^ I 2 ' 



(27) 



In the thermodynamic limit, the key constants Fi and F2 
with (S", s) = (1, 1/2) are estimated to be 0.304887 and 
-0.337779, respectively. Up to order 0(5""), we end up 
with the Hamiltonian 



where 



Eg — Ec\a 



Eo + Ei 



(28) 



(29) 
(30) 



In Fig. ^ we plot usj^ and cuf as functions of k at 
(S,s) = (1,1/2) with the previous numerical estimates 
||l4| for finite chains. We find that the antiferromagnetic 
mode is now improved to a great extent. Furthermore 
the band widths of the two branches, which are exactly 
the same within the noninteracting spin- wave theory, are 
now indeed different, due to the interactions. The gap 
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A = J is replaced by Z\ = (1 - 2^2) J ~ 1.676 J, which is 
much closer to the exact value 1.759J. The ground-state 
energy Eg is also refined: i^ciass + Eq ~ —1.437 J, while 
-E-ciass + Eq + El ~ —1.459 J, where the exact value is 
-1.454J. 



B. Thermodynamics 

Now let us start our MSW theory from the above- 
obtained improved dispersion relations. At finite tem- 
peratures we replace a^afc and PlPk m the spin-wave 
Hamiltonian by = J2n- n+ ''^^-Pk(ji~ ,n~^), where 
Pk{n^ ,n^) is the probability of n~ ferromagnetic and 
nf^ antiferromagnetic spin waves appearing in the k- 
momentum state and satisfies 



^ Ffe(n-,n+) = l 



(31) 



for all fc's. Then the free energy at zero field is given by 

F = E,. 



fceT^ ^ Pk{n ,n+)lnPfc(n , 



(32) 



We now carry out the minimization of the free energy 
( ^ ) with respect to Pfc(n~, n"'")'s under a particular con- 
straint as well as the trivial constraint (p^) . The original 
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idea introduced by Takahashi was that the zero- 

field magnetization should be zero. This constraint works 
quite well especially in ferromagnets, where it serves to 
control the number of Holstcin-PrimakofF bossons. Let 
us apply the same constraint to the present model: 



= (5-,s)7V-^^anr =0, (33) 



k <T=± 



where — J2j Sj and = J^j *j ■ Equation ( p3D claims 
that the thermal fluctuation <^'^k'^ should be con- 

strained to take the classical magnetization (5 — s)N. 
The free energy and the magnetic susceptibility at ther- 
mal equilibrium are then given by 

F^Eg+fi{S^ s)N - fceT^ ^ ln(l + n^) , (34) 



k u=± 



^^n-(l-fn-), 



fe t7=± 



with 



1 



(35) 



(36) 



where fi is a Lagrange multiplier due to the condition 
(^). The susceptibility has been obtained by calculat- 
ing X = igfiB)H{M^) - (M)2)/3iVr M. Equations (H) 
and (^) are expanded in powers of at low temper- 
atures as 



32{S - s)2{Ssy- 



{S^ + Ss + s^)al) 4C(i) 



+0{t 



72\ 



-t"" - iSs)^{S - s)5 2 + (5 - s) 



2tt 



ah: 



<-i + 0(t-5) 



(37) 
(38) 



where C{z) is Riemann's zeta function and t = k^T / J = 
fceT/JT with7 = \ — ri{S-\-s)/ Ss — r2/ Ss. Surprisingly, 
the low-temperature series expansions (^7|) and (|3|) are 
exactly the same as the thermodynamic Bethc-ansatz cal- 
culations for the spin- 1/2 ferromagnet js^ except for 7. 
Thus we recognize similarities between ferrimagnets and 
ferromagnets at low temperatures. We note, however, 
that ferrimagnets of 5* = 2s should not strictly be identi- 
fied with spin-s ferromagnets because of the scaling factor 
7. With the interactions, the original temperature T is 
replaced by T(> T) or equivalently the original spins are 
reduced. Therefore ferrimagnets of 5 = 2s behaves like 
spin-s^(< s) ferromagnets, at least at low temperatures, 
which reminds us of the quantum spin reduction jl2j in 
ferrimagnets. The spin- wave theory shows that the stag- 
gered magnetization is reduced to {S + s)N — 2t with 
2T = Y.^[{S + s)/suk-l]. 

Unfortunately, in the present model, the constraint in- 
troduced above is not useful at all at high temperatures 



because it allows the number of bosons of each mode to 
diverge. Actually, under the condition of zero magnetiza- 
tion, the specific heat becomes a monotonically increas- 
ing function. Hence we propose an alternative constraint 
|17| . Let us consider the minimization of the free energy 
constraining the staggered magnetization to be zero: 

(: 5- - S-) = (5 + s)7V - (5 + s) ^ ^ ^ ^ , 

k a=± ^'^ 

(39) 

where the normal ordering is taken with respect to 
both operators a and (3. Equation ( |3^ ) implies, just 
as Eq. ( ^3|) does, that the thermal fluctuation [S + 
s) nf,/ujk should take the c/assica/ value (S+s)N^ 

rather than the renormalized value {S + s)N — 2t. We 
note that, based on the naive idea of zero staggered mag- 
netization without the normal ordering, we have 
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{S + s)N^2t-{S + 



(40) 



Besides such a phenomenological argument, a more strin- 
gent reason drives us to adopt Eq. (p9|): The condition 
( |39| ) exactly reproduces the low-temperature series ex- 
pansions ( p7|) and (|38|), whereas the condition ( ^o| ) fails 
to do. The conventional spin-wave theory at least brings 
the correct leading term of the specific heat, which is 
obviously reproduced by the MSW theory with the con- 
straint of zero magnetization. The low-temperature be- 
havior must be inherent in the dispersion relation and 
should not be influenced by supplementary constraints. 
Now the set of self-consistent equations (^) and ( |35| ) is 
replaced by 



F^E^ + fiiS + s)N -kBTY^Yl 1^(1 + ^fe) ' (41) 

(42) 



k <y=± 



3A;bT 



fe <T=± 



with 



1 



^lJuj^~-tJ.(S+s}/<^k]/kBT _ 2 



(43) 



where fi is a Lagrange multiplier due to the condition 

In the case of (S, s) — (1,1/2), we have numerically 
obtained ( pl| ) and (^3) in the thermodynamic limit and 
visualize them in FigT^ and Fig. ^ where the solid curves 
represent the calculations starting from the improved dis- 
persion relations ujf , while the dashed curves the calcula- 
tions with LU^ [|l^ instead of . Though the present cal- 
culation including the 0(5°) interactions underestimates 
the height of the Schottky-like peak, the interactions in- 
deed correct the location of the peak, which emphasizes 
that the gapped antiferromagnetic spin-wave mode pre- 
dominantly contributes to this peak. We would like to at- 
tribute the discrepancy unsolved to the constraint rather 
than to the dispersion relations. Even though we adopt 
the two constraints ( |33| ) and (^) simultaneously, the re- 
sults do not change essentially. This is not so surprising 
because the two constraints play almost the same role 
at low temperatures, as the low-temperature series ex- 
pansions imply, whereas only the constraint ( |39| ) is rel- 
evant at high temperatures. The precise description of 
the overall temperature dependence might be obtained 
with a temperature-dependent constraint, which is not 
so interesting. 

We obtain the best results from the MSW theory at low 
temperatures. Figure ^(b) fully convinces us of the valid- 
ity of the present MSW calculation. Figure ||(b) clearly 
reveals the low-temperature behavior which no numeri- 
cal tool has succeeded in observing. Takahashi compared 
his MSW findings for ferromagnets with the spin- 
1/2 thermodynamic Bethe-ansatz calculations and 



found that the MSW theory correctly describes the lead- 
ing two terms of the specific heat and the leading three 
terms of the susceptibility in its low-temperature series 
expansions. We also find that the low-temperature series 
expansions (^) and (^8|) coincide with ones of the spin- 
1/2 ferromagnet to the same extent except for the scal- 
ing factor 7. That may be why the MSW and numerical 
findings for the susceptibility join at higher temperatures 
than those for the specific heat. 



IV. SUMMARY 



Developing an analytic argument as well as employing 
various numerical tools, we have investigated the ther- 
modynamic properties of Heisenberg ferrimagnetic spin 
chains. Both ferromagnetic and antiferromagnetic as- 
pects lie in the model; they are most clearly exhibited 
at low and intermediate temperatures, respectively. One 
might say that mixed-spin chains possess mixed features. 
We have shown that the MSW theory starting from the 
improved dispersion relations precisely describes the low- 
temperature behavior of the model. We appeal to ex- 
perimentalists to carry out specific heat and susceptibil- 
ity measurements of mixed-spin materials especially at 
low temperatures. On the other hand, the Schottky- 
like peak of the specific heat and the minimum of the 
susceptibility-temperature product clearly reflecting the 
antiferromagnetic gap imply that the antiferromagnetic 
excitations are not smeared out in the ferromagnetic 
spectra but stand out clearly. Hence neutron-scattering 
measurements are also encouraged. 
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FIG. 1. Temperature dependence of the specific heat per 
unit cell: (a) QMC findings (Q) and DMRG results (x) at 
oo. The dotted line represents the Schottky-type spe- 
cific heat (2.9). The numerical uncertainty is smaller than 
the symbol size, (b) QTM calculations at various values of 
A^, where we show the A — > cxa curve as well, which is ob- 
tained by interconnecting the QMC findings at ksT/J > 0.4 
and the DMRG results at ksT/J < 0.4. 



FIG. 2. Schematic representations of the A/ = N/2 ground 
state of the A = 1 isolated dimer composed of spin 1 and 
spin 1/2 and its ferromagnetic (b) and antiferromagnetic (c) 
excitations. The arrow (the bullet symbol) denotes a spin 1 /2 
with its fixed (unfixed) projection value. The solid (broken) 
segment is a singlet (triplet) pair. The circle represents an 
operation of constructing a spin 1 by symmetrizing the two 
spin 1/2 's inside. 



FIG. 3. Temperature dependence of the magnetic suscepti- 
bility times temperature per unit cell: QMC findings (Q) and 
DMRG results (x) at N oo. The numerical uncertainty is 
smaller than the symbol size. 



FIG. 4. Dispersion relations of the lowest-energy states in 
the subspaces of M = 1: The noninteracting spin- wave 

result (dotted lines) and an improved calculation taking into 
account interactions between spin waves (solid lines). Previ- 
ous numerical calculations [14] are also shown for the sake of 
comparison. Here we plot the excitation energy E{k) taking 
the ground-state energy and twice the lattice constant as zero 
and unity, respectively. 



FIG. 5. Temperature dependence of the specific heat per 
unit cell: The MSW calculation with noninteracting spin 
waves (dotted lines) and that with the improved dispersion 
relations (solid lines). Numerical findings [QMC (O) and 
DMRG(x)l are also shown. 
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FIG. 6. Temperature dependence of the magnetic suscep- 
tibility times temperature per unit cell: The MSW calcula- 
tion with noninteracting spin waves (dotted lines) and that 
with the improved dispersion relations (solid lines). Numeri- 
cal findings [QMC (O) and DMRG(x)] are also shown. 
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